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Abstract 

Using the Fokker-Planck equation we have studied the drag co-efficient A(t) and the conse- 
quent shift Apj_(L) in the transverse momentum due to collisional energy loss of energetic 
partons while passing through a chemically equilibrating quark-gluon plasma. Based on 
these we estimate the quenching factor Q(p±) when the medium is undergoing longitudi- 
nal expansion governed by master rate equations. In contrast to the case of chemically 
equilibrated plasma investigated earlier by Mustafa and Thoma [2] we find less quenching 
because our calculated Q(p±) is always greater at all momenta. This result is attributed 
to the weak drag coefficient operating during initial state interactions. 



PACS: 12.38.Bx; 24.85.+p 
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1 Introduction 



The inclusive yield of hadrons produced with high transverse momentum p± in Au+Au 
collisions at the Relativistic Heavy Ion Collider (RHIC) has recently been shown to be 
significantly suppressed in comparison with the cumulative yield of binary NN interactions. 
This effect, called "jet quenching", is considered an important signal for the production of 
a quark-gluon plasma (QGP) and a vast literature [1]-[21] already exists dealing with the 
phenomenological/theoretical aspects of the said effect; see Appendix for a brief review. 

Jet quenching is believed to occur due to the fact that hard partons produced at the 
initial stage of the heavy ion collision lose energy as they propagate through the fireball 
making Q(p±) < 1. Here the quenching factor is experimentally defined by the ratio 

where is the inclusive hadron spectrum at transverse momentum p±, and the su- 

perscripts refer to the medium and vacuum, respectively. The theoretical calculation of 
Q invokes two important, mutually competing, mechanisms of the energy loss described 
below. 

The radiative mechanism [21] is caused by deceleration of the colour charge accompanied 
by the bremsstrahlung of soft gluons; perturbative QCD allows the determination of the 
soft gluon distribution I{uj) in terms of a characteristic frequency u c . 

On the other hand, the collisional mechanism [9]- [17] arises from the elastic encounters 
with the other partons of the medium; here a QCD motivated drag force permits the 
evaluation of the average momentum loss Ap± suffered by the test parton over a specified 
distance L. The collisional loss theory has been recently applied by Mustafa and Thoma [2] 
(referred to as MT hereafter) assuming a temperature cooling law relevant to chemically 
equilibrated QGP. The aim of the present paper is to extend their theme to the case where 
the evolving fugacities have not yet achieved chemical equilibrium. 

For the sake of convenience Sec. 2 below defines our notations and briefly recapitulates 
the main formulae derived by MT. The details of our numerical work are presented in Sec. 3. 
Finally, physical interpretations of the results along with some concluding remarks appear 
in Sec. 4. 



2 Recapitulation of MT Theory [2] 

Step i) In the fireball rest frame a Taylor expansion of the hadron spectrum £ is made [21] 
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to rewrite (1) in the standard form 

S med (p ± ) w S vac (p ± + AE); AE = J dee D{e) (2) 

where e is the random collisional energy loss over time span t, E = E + e the random 
surviving energy, D(e) the probability distribution that a parton loses the energy e, and 
AS the average energy loss in traversing the distance L = c t. 

Step ii) Let the instantaneous momentum of the leading parton moving in the transverse 
x direction be called p —\ p | e x . Assuming spatially uniform plasma a Boltzmann transport 
equation is set-up for the distribution function D(t,p), Landau's approximation is made in 
the collision integral, and the resulting Fokker-Planck equation is cast in the form 

3D 8 r ->, _i d 2 



dt 



T x ^p)D + — [T 2 (t,p)D]. (3) 



—* 

The transport coefficients T\ and T 2 are defined in terms of the net instantaneous collision 
rate w(t,p, k), involving soft momentum transfer k, by 

fi(t,p) = f d 3 k w(t,p,k)k (4) 



r 2 (t,p) = -Jd"kw(t,p,k)k 2 (5) 

w(t,p,k) = lij J^sffi'fi Vicl a i> ^ 



where j labels different species (quarks, antiquarks, gluons) present in the plasma, 7j is the 
degeneracy factor, fj the species' distribution function depending implicitly on the time 
through the fugacity Xj(t) and temperature T(t), v TC \ the relative speed between the test 
parton and j-th species, and Uj the associated scattering cross section. 



Step Hi) Correspondence with classical dissipative motion tells that the function T\ 
scales like p and represents the rate of energy loss —dE/dx. Hence, a mean drag coefficient 
A(t) is constructed in one-dimensional notation via 



1 dE 



A(t) = (-- S = / d 3 p =1M= / [ d 3 p ^J= (7) 

p dV J expy/p 2 + m2/T- 1 J exp ^jp 2 + m 2 g /T-l 

Since the QGP expected at RHIC and LHC is likely to be out of chemical equilibrium 
it is necessary to investigate the energy loss in this case [23]. Indeed, even away from 
chemical equilibrium, dynamical screening remains operational within the HTL-resummed 
perturbation theory. More explicitly, the collisional energy loss for a heavy quark (mass M) 
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propagating through a QGP parametrized in terms of the distribution functions \ q np and 
\ g riB, respectively, where \ q>g are the fugacity factors describing chemical non-equilibrium, 
becomes [1] 



dE . 2 (1 + 9/4) 



0.920^?^ 2 A " JV ' /(12A9+2A * JV ' ) 



(8) 



This expression [1, 24] is valid for energetic quarks with E 3> M 2 /T and contains for 
\ q = \ g = 1 the original result of [2]. The screening mass parameter is 



mi = ATia s {\ g + X q N f /Q)T 2 /3. 



(9) 



The expression reduces to the expression given in Mustafa and Thoma paper by putting 

-9/ 
2 



1. It should be noted that the factor ( 1+ y 4 ) arises due to averaging over the 



quark and gluon contributions. 

Next, the analogy with Einstein's random walk relation is exploited to identify the 
momentum-averaged diffusion coefficient 



V F (t) = (T 2 (t,p))=AT 2 



(10) 



so that the evolution equation for the momentum-distribution of the Brownian particle 
becomes 

dD A d . . d 2 D 

- = A - m+Vrw (11) 

Step iv) Employing Fourier transform and method of characteristics the above partial 
differential equation is solved analytically subject to the initial condition D(p, 0) = 5{p— p ). 
The result in terms of length and energy variable is 

I ! E — En ■ B(L) X 2 

D(L, E) 



cxp 



(E-E -B(L) 



W(L) 



(12) 



where 



B(L) = exp ^ - ^ dt' A{t') j , 



and 



W(L) = 4 / rft' V F (t') exp 2 / di" A(t") B 2 (L) . 
Jo Jo 

Step v) The mean energy (E) and average energy loss AE are determined from 

(E) = (p) = / dE E D(L, E) = E B{L) 
Jo 
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(13) 



(14) 



(15) 



and 

AE = Ap = E - ( E ) = E ( 1 - B(L) ) (16) 

Eqs.(13, 14) can now be inserted into the basic expression (2) along with the following 
high-energy parametrization for jet hadronization at RHIC: 

£ Va V)= const (l + P fJ V , (17) 
with v = 8 and P n = 1.75 GeV. 



Step vi) Finally, the geometry of the head-on heavy ion collision is accounted for by 
considering a cylindrical plasma of radius R and the test particle moving in the central 
rapidity region. If the latter was created at the location (r, 0) in the transverse plane z = 
then it travels a distance 

L(r,<f>) = (R 2 - r 2 sin 2 0) 1/2 - r cos0 (18) 

before leaving the cylinder. Upon averaging (2) over the creation configuration the effective 
quenching factor becomes 

Q(P±) = fj d<f> J* d 2 r S VaC (p ± + Ap) I S VaC (p ± ) (19) 

which is amenable to direct computations. 



3 Numerical Application 

Mustafa- Thoma procedure: MT assumed the QGP to expand longitudinally according to 
Bjorken's boost-invariant hydrodynamics [26] so that the temperature T(t) on the trans- 
verse plane decreases according to the scaling law 

T(t) = T(t )(^) 1/3 , (20) 

where t the instant when the background plasma had just attained local kinetic as well 
as chemical equilibrium so that all partonic fugacities had become unity. For an A + B 
collision at RHIC, MT took 

to = 0.3 fm; T = 0.5 GeV (21) 
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Table 1: Different sets of initial conditions of the temperature, fugacities and parton number 
densities at r = 0.7fm/c for RHIC and r = 0.5fm/c for LHC. 

RHIC(l) LHC(l) RHIC(2) LHC(2) RHIC(3) LHC(3) 



T(GeV) 


0.55 


0.82 


0.55 


0.82 


0.4 


0.72 


A, 


0.05 


0.124 


0.2 


0.496 


0.53 


0.761 


\ 


0.008 


0.02 


0.032 


0.08 


0.083 


0.118 


n g (fm 3) 


2.15 


18 


8.6 


72 


8.6 


72 


n g (fm~3) 


0.19 


1.573 


0.76 


6.29 


0.76 


6.29 



The momentum averaging in (7) was done using a Boltzmann distribution at temperature 
T(t) for gluons. 

Our Procedure: For the case of longitudinal expansion we know that, in the early stage 
of evolution, the plasma may achieve equilibrium thermally but not yet chemically. Par- 
tonic reactions drive the system towards chemical equilibration through the master rate 
equations [28] 

£ + 3| + i = «,(l-A.)-2ft(l-3), 

i + 3 i + I = W^-f), 

A, T r Oi \A ? Aj/ 

/ b \ 3/4 

|A 9 + ^AJ T 3 r = const, (22) 

where the usual meaning of the symbols can be found in Ref. [28]. Their solutions subject 
to the initial conditions [29] (Table 1) give the temperature and fugacities at general t, 
and these enter the distributions fj(t,q) of various species written in (6). But the mo- 
mentum integration in (7) was performed employing an equilibrium Bose-Einstein form at 
temperature T(t). 

Computed Results: Fig.l plots the collisional drag coefficient A versus elapsed time t for 
energetic gluons, in the QGP phase of the expanding fireball, obtained from (7). The 
dashed-dotted line corresponds to MT's equilibrated plasma described by the cooling law 
(20) whereas other curves refer to our equilibrating plasma governed by the master equa- 
tions (22). Next, Fig.2 shows the relative energy loss Ap±/p±_ versus the transverse distance 
L. Finally, Fig.3 displays the effective quenching factor Q versus the transverse momentum 
p±. Logical explanation of these results is taken- up in the next section. 



6 




Figure 1: The drag coefficient A(t) in an expanding QGP. In the left panel, the solid, 
dotted and dashed line denotes chemically equilibrating QGP corresponding to RHIC(l), 
RHIC(2) and RHIC(3) initial conditions, respectively. The dashed-dotted line corresponds 
to MT's equilibrated system. On the other hand right panel corresponds to LHC energy. 
MT choose the initial values from eq.(21) with X q = X g = 1 whereas latter takes the initial 
conditions from Table 1 [29]. 

4 Interpretations and conclusions 

As is well known the distribution function fj(t, q) of the j-th species in (6) grows with higher 
temperature T(t) and fugacity Xj(t). The numerical solution of (20) and (22), with the 
initial conditions stated, reveals that Mustafa-Thoma fugacities (and hence the transition 

— # 

rates w(t,p, k)) are about an order of magnitude more than those in our work at all times 
of experimental interest. Consequently, in the left panel of Fig.l the drag coefficient A MT 
is significantly greater than A RHIC(1) ; indeed the ratio A MT /A Kmc( -^ ps 2 at t « 1 fm. The 
reason for this difference lies in the fact that MT deals with a fully equilibrated system 
whereas RHIC(l) deals with the highly unequilibrated system as is evidenced from their 
initial conditions. As a consequence, gluon density in MT system becomes larger than 
RHIC(l) resulting higher value of the drag coefficient A(i). At large time (greater than 
8 fm), our result does not differ much from MT result because as time elapses chemically 
unequilibrated system (RHIC(l)) approaches towards chemical equilibrium making the two 
system alike. However, our calculated value is always less than A MT . Thus equilibrating 
QGP provides less drag force in comparison to equilibrated QGP. 

However for the illustration of the dependence of A(t) on initial parton fugacities, 
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Figure 2: The effective shift of the scaled transverse momentum Ap±/p± as a function of 
traversed distance L. The notations of the curves are the same as in Fig.l 

we have calculated A(t) at RHIC(2) and RHIC(3) [29] where the values of A(t) is not 
much different from MT. This is due to the fact RHIC(2,3) is not much unequilibrated as 
RHIC(l) as is evidenced by their initial gluon/quark fugacities or gluon densities compared 
to RHIC(l). 

In the right panel of Fig.l the drag coefficient A(t) at LHC energy also decreases with time. 
But our result at LHC energy is always somewhat greater than the value at RHIC. Thus 
the drag coefficient increases with increase in the center-of-mass energy. 

Next, in the left panel of Fig. 2 the relative energy loss Ap/p ± | MT far exceeds Ap/p ± |° ur 
computed at sizable distances; indeed their mutual ratio is about 1.4 at L ps 7 fm at 
RHIC(l) energy. This happens because for higher integrand A(t') the function B(L) be- 
comes lower in (13, 16). Of course, in both cases Ap scales almost linearly with p± for a 
given L. In the right panel of the Fig.2 the relative energy loss Ap/p± |° ur at LHC energy 
does not differ from the RHIC value. Thus increase in center-of-mass energy in going from 
RHIC to LHC the relative energy loss in equilibrating QGP will not be affected. 

Finally, in the left panel of Fig. 3, the quenching factor Q(p±) is stronger at lower p± in 
our calculation as well as MT calculation and gradually weakens at higher p± in agreement 
with experimental results [20]. However, in our case, the quenching factor Q(p±) calculated 
for equilibrating QGP shows less suppression in comparison to quenching factor Q(p±) 
calculated for equilibrated QGP. This is again a reflection of the fact that, in contrast to 
MT equilibrated case, our energy loss Api± is much smaller over the whole assembly of 
lengths L encountered in (18, 19). Moreover, the qualitative behavoiur is same in both the 
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Figure 3: The quenching factor Q(p±) as a function of transverse momentum p±. The 
notations of the curves are the same as in Fig.l 

cases. In right panel of the Fig. 3, the quenching factor Q(p±) shows stronger suppression 
at LHC energy in comparison to RHIC energy at lower values of the transverse momentum 
less than p± 8 GeV. 

We end the paper with two concluding remarks. Although our assumptions of chemi- 
cally equilibrating QGP is more realistic than MT's assumption of equilibrated plasma yet 
the equilibrating QGP will provide weak drag force for the jets in comparison to equili- 
brated QGP. Thus most jet suppression occurs after the equilibration of the QGP has been 
achieved. However, this initial state effect cannot be neglected. Therefore, jet quenching 
cannot be exactly called as a final state effect. 



ACKNOWLEDGEMENTS 

BKP is thankful to V.J. Menon for the useful discussion during this work. BKP is also 
thankful to B. Miiller, M. Thoma and M.G.Mustafa for their useful private communication. 
BKP acknowledges the financial support to IIT Roorkee for the Faculty Initiation Grant 
and MM thanks to CSIR, New Delhi for his financial support. 



9 



References 

[1] R. Baier, D. Schiff, B.G. Zakharovo, Ann. Rev. Nucl. Part. Sci. 50 (2000) 37. 

[2] M. G. Mustafa and Markus Thoma, Acta Phys. Hung. A 22 (2005) 93. 

[3] M. G. Mustafa, Phys. Rev. C 72 (2005) 014905. 

[4] M. Gyulassy and M. Pliimer, Phys. Lett. B 243 (1990) 432. 

[5] J.D. Bjorken, Fermilab preprint 82/59-THY (1982, unpublished). 

[6] B. Svetitsky, Phys. Rev. D 37 (1988) 2484. 

[7] S. Mrowczyhski, Phys. Lett. B 269 (1991) 383. 

[8] Y. Koike and T. Matsui, Phys. Rev. D 45 (1992) 3237. 

[9] E. Braaten and R.D. Pisarski, Nucl. Phys. B 337 (1990) 569. 
[10] E. Braaten and M.H. Thoma, Phys. Rev. D 44 (1991) 1298. 
[11] E. Braaten and M.H. Thoma, Phys. Rev. D 44 (1991) 2625(R). 
[12] H. Vija and M.H. Thoma, Phys. Lett. B 342 (1995) 212. 
[13] P. Romatschke and M. Strickland, hep-ph/0309093. 
[14] M.H. Thoma and M. Gyulassy, Nucl. Phys. B 351 (1991) 491. 
[15] M. H. Thoma, Phys. Lett. B 273 (1991) 128. 
[16] M. H. Thoma, J. Phys. G 26 (2000) 1507. 
[17] M. H. Thoma, Eur. Phys. J. C 16 (2000) 513. 

[18] R. Baier, D. Schiff, and B.G. Zakharov, Ann. Rev. Nucl. Part. Sci. 50 (2000) 37. 
[19] E. Wang and X.-N. Wang, Phys. Rev. Lett. 87 (2001) 142301. 

[20] K. Adcox et al. [PHENIX Collaboration], Phys. Rev. Lett. 88 (2002) 022301; C. Adler 
et al. [STAR Collaboration], Phys. Rev. Lett. 89 (2002) 202301; B.B. Back et al. 
[PHOBOS Collaboration], nud-ex/ 0302015. 

[21] R. Baier, Y. L. Dokshitzer, A. H. Mueller, and D. Schiff, JHEP 0109 (2001) 033. 

[22] B. Miiller, Phys. Rev. C 67 (2003) 061901(R). 

10 



[23] Mustafa MG, Pal D, Srivastava DK., Phys. Rev. C57 (1998) 889. 

[24] Baier R, Dirks M, Redlich K., hep-ph/9910353 (1999). 

[25] Braaten E, Thoma MH., Phys. Rev.D44 (1991) R2625. 

[26] J. D. Bjorken, Phys. Rev. D 27 (1983) 140. 

[27] R. J. Fries, B. Midler and D. K. Srivastava, Phys. Rev. Lett. 90 (2003) 132301. 

[28] T. S. Biro, E. van Doom, M. H. Thoma, B. Midler, and X.-N. Wang, Phys. Rev. C 48 
(1993) 1275. 

[29] X.-N Wang and M. Gyulassy, Phys. Rev. D 44, 3501 (1991). 



11 



Appendix: Jet quenching summarized 



The yield of hadron produced with high transverse momentum pj_ in Au+Au collisions at 
RHIC has recently shown to be significantly suppressed in comparison with the cumulative 
yield of NN collisions. This effect, so called "Jet quenching" was predicted to occur due to 
energy loss suffered by hard scattered partons. The energy loss is expected to occur due to 
interaction of the hard partons with the surrounding dense medium. 

There are two contributions to the energy loss of the patons in the medium. One is due 
to the collisions among the partons in the medium and other due to the radiation emitted 
by the decelerated colour charge i.e., bremsstrahlung of gluons. 

In the initial stage of ultrarelativistic collisions energetic partons are produced from 
hard collisions energetic partons are produced from hard collisions between the partons of 
the nuclei. Receiving a large transverse momentum, these partons will propagate through 
the fireball which might consists of quark-gluon plasma phase for a transitional period of 
about few fm/c. These high energy partons will manifest themselves as jets leaving the 
fireball. These energy partons will loose energy due to interaction of the hard partons with 
the fireball medium. Hence jet quenching will result. The amount of jet quenching might 
depend on the state of the fireball i.e., QGP or hot hadron gas, respectively. Therefore, jet 
quenching has been proposed as the possible signature of the QGP formation. 

It is difficult to measure the energy loss of the scattered partons directly in heavy-ion 
collision because of the large multiplicity of the emitted hadrons makes it almost impossible 
to isolate the resulting jet by the kinematic cut. However, this energy loss of the hard 
partons may affect the equivalent loss of the energy of the leading hadrons produced it its 
fragmentation, this is what has been observed at RHIC experiment []. Preliminary data 
from run 2 at RHIC Au+Au at ^snn = 200 GeV confirm the effect observed in run 1 and 
its interpretation as jet quenching. For pions with p±_ pa 5 GeV/c the measured suppression 
factor is about 1/5. 
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